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Abstract 

We report an improved measurement of the branching fraction for D*j(2317)^if and 

present evidence of the Dsj{2AQQ)'^ K~ decay. These results are obtained from a data 

sample containing 386 million BB pairs that was collected near the T(4S') resonance, with the 
Belle detector at the KEKB asymmetric energy e+e" collider. 

PACS numbers: 13.25.Hw, 14.40.Nd 



3 



Two narrow resonances denoted as D*7{'2,317)~^ and Dsj{2460)~^ have been observed re- 
cently in e'^e~ continuum interactions ^ decays 0|- The surprising! 



low masses and small widths of these states initiated a wide theoretical discussion 
Although the O"*" and l"*" quantum numbers have been established for the D*j{2317)~^ and 
-Dsj(2460)'^ resonances 0, respectively, the nature of these states is still unclear. 

In this paper we report an updated study of the decays 5° D^jK" with a data 
sample that is approximately 2.5 times larger than in the recently Belle published pa- 
per that first reported the — > D*j(2317)~^i^r~ decay mode. In the previous Belle 
study the product branching fraction B(B^ D*j {2317)+ K-)x 13 (JJ*j (2317)+ D+n^) = 
{5.3t\i ± 0.7 ± 1.4) X 10"^ was measured and an upper limit B{B° Dsji2460)+K ) x 
-B(D^j(2460)+ < 0.94 x 10"^ was set. These measurements [ld\ show that 

B{W D*j{2317yK-) is of the same order of magnitude as ^ D+R-) [U E3 

and at least a factor of two larger than the branching fraction for D sj{2AQ{))+ K~ . 

The B^ — »• D+^j^K~ decays can be described by a PQCD factorization W exchange 
process or, alternatively, by final state interactions [3, [l^ . Assuming there is a 

four-quark component of the Dgj mesons, the tree diagram with ss pair creation may also 
contribute Although accurate theoretical calculations of branching fractions are difficult 
for these decay modes, the experimental results disagree with the naive expectation that 
the ratio B{W -* D+h~)/B{W Djjh-) should be similar for h' = vr", K' or D^- 

This analysis is based on a large data sample, which contains 386 million BB pairs, 
collected with the Belle detector at the KEKB asymmetric-energy e+e~ (3.5 on 8 GeV) 
collider jl^ operating at the T(45') resonance. The Belle detector is a large-solid- angle 
magnetic spectrometer that consists of a silicon vertex detector (SVD), a 50-layer central 
drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like 
arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) located inside a super-conducting solenoid coil that 
provides a 1.5 T magnetic field. An iron fiux-return located outside of the coil is instru- 
mented to detect mesons and to identify muons (KLM). The detector is described in 
detail elsewhere [Til]. Two inner detector configurations were used. A 2.0 cm beampipe and 
a 3-layer silicon vertex detector was used for the first sample of 152 million BB pairs, while 
a 1.5 cm beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were used 
to record the remaining 234 million BB pairs |2ol|. 

In this analysis we applied the same selection criteria as in 0], where a detailed description 
of the criteria can be found. The only differences between the two analyses arise due to the 
vertex detector upgrade. According to the MC simulation, minor differences in the signal 
widths and efficiencies are expected for the two SVD subdetector configurations, leading to 
respective corrections applied in the fit procedure and efficiency calculations. 

Kaon and pion mass hypotheses are assigned to the charged tracks with momenta p > 
lOOMeV/c llOj using a likelihood ratio Cr/tt = '^k/{'^k + '^■k), obtained by combining 
information from the CDC (dE/dx), ACC, and TOF systems. We require Cx/n > 0.6 
{Ck/w < 0.6) for kaon (pion) candidates [19^ . 

ECL clusters with a photon-like shape and energies larger than 50 MeV, that are not 
associated with charged tracks, are accepted as photon candidates. Photon pairs of invariant 
mass within ±12MeV/c2 (~ 3a in the 7r° mass resolution) of the 7r° mass are considered 7r° 
candidates; the 7r° momentum is required to be larger than lOOMeV/c. 

Kg candidates are formed from tt+tt^ pairs with an invariant mass within ±10MeV/c^ 
(~ 3a) of the nominal Kg mass. Invariant masses of K*^ K+tt~ candidates are required 
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to be within ±50 MeV/c^ of the nominal K*^ mass; those of — 5> K^K~ candidates, within 
±12MeV/c^ of the mass. mesons are reconstructed in the 0vr+, K*^K^ and KgK^ 
decay channels; a mass window of ±12MeV/c^ (~ 2.5cr) is imposed in each case. The 
Dgj mesons are reconstructed in the D*j(2317)+ -D^vr° and Dsj(2460)+ D'^l decay 
modes; within the mass difference ranges \M(Dfn^) — M(Df) — 348. 6| < 20MeV/c^ and 
|M(D+7) - M{p+) - 487.91 < SOMeV/c^. 

Candidate — > DfjK~ and D^jTc^ are formed and the signal is extracted us- 
ing the energy difference AE = i?^*^ — -Ebeam beam-constrained mass Mbc = 



(-^beam)^ ~ (Ps^^)^) -E^*^ and are the energy and momentum of the B candidate 
in the center-of-mass (CM) system and -Ebeam beam energy. Only events within 

the intervals Mbc > 5.2 GeV/c^ and \ AE\ < 0.2 GeV are used in this analysis. The B meson 
signal region is defined by \AE\ < 0.04 GeV and 5.272 GeV/c^ < Mbc < 5.288 GeV/c^. 

Combinatorial background for channels involving the Dsj(2460)"'" was further suppressed 
by requiring cos6'£)^-y < 0.7. The hehcity angle Ods-t defined as the angle between the di- 
rection opposite the B momentum and the momentum in the -D^7 rest frame. This re- 
quirement rejects 49% of background events and only 6% of signal events, assuming 
for the Dsj(2460)+. The uncertainty due to this assumption is included in the systematic 
error. 

For events with two or more B candidates, the and vr" candidates with invariant 
masses closest to their nominal values and the B daughter or vr" candidate with the 
best Ck/tt value are chosen. No multiple entries are found in the data. 

We exploit the event topology to separate BB events (spherical) from the continuum 
background (jetlike). The ratio of the second and zeroth Fox- Wolfram moments jlJ] of all 
particles in the event is required to be less than 0.5. For such events, we form a Fisher 
discriminant from six modified Fox- Wolfram moments. A signal (background) likelihood £s 
[Cbg] is obtained using signal MC (sideband) data from the product of probability density 
functions for the Fisher discriminant and cos6b, where 6b is the B flight direction in the CM 
system with respect to the z axis. We require 71 = Cs/{Cs + Cbg) > 0.4 for K*^K^ 
and TZ > 0.25 for the other decay modes, which have lower backgrounds. 

The AE and AM{Dsj) distributions for the DfjK~ combinations are shown in Fig. ^for 
the range 5.272 GeV/c^ < Mbc < 5.288 GeV/c^. To obtain the AM^D^j) distributions we 
relax the AM{Dsj) requirements and apply the tight selection on AE. The AE distributions 
are modelled using a linear background function and a Gaussian signal shape (the Crystal 
Ball shape function i^] is used for the Dgj (2460)^) with zero mean and a fixed width 



determined from MC data. The AM{Dsj) distributions are described by the sum of a 
signal Gaussian and a linear background. The widths of the Gaussians are fixed from MC 
while the peak positions are allowed to float. A strong B^ — > D*j{2317)^K~ signal is 
observed and evidence of the B^ — Dsj{2460)^K^ signal is also seen. The Gaussian peak 
positions obtained from the fits correspond to the Dgj mass values of 2319.2 ± 1.3MeV/c^ 
and 2456.2 ±6.5 MeV/c^ for the D*j(2317)+ and D,j(2460)+, respectively. These values are 
in a good agreement with the most recent BaBar measurements ji^l of Dgj masses in the 
continuum, 2318.9 ± 0.3 ± 0.9 MeV/c^ and 2459.4 ± 0.3 ± 1.0 MeV/c^. 

Signal yields, efficiencies, branching fractions and significances for the studied decay 
channels are shown in Table 1. The signal yields are obtained from the fits of histograms 
shown in Fig. Q where the three Dg decay channels are combined. The B^ — > D^jK~ 
branching fractions and significances are obtained using a simultaneous fit to the AM{Dsj) 
distributions for the three decay channels, with independent background descriptions. 
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FIG. 1: AE (a) and AM{Dsj) (b) distributions for the 5° ^ D*j{23l7)+K~ decay, and AE (c) 
and AM{Dsj) (d) distributions for the W Dsj{2Am)+K- decay. 



but common values for the signal width (fixed from MC) and peak position (allowed to float). 
The branching fractions obtained in the individual decay modes agree within the statistical 
errors. The three error terms are the statistical uncertainty, the total systematic error, and 
the uncertainty due to branching fractions; this last term is dominated by the ~ 25% 
uncertainty in B{D^ (pn^) [24,J. The last two systematic terms are combined for the 

5° Dsj{2AQ0)^K~' decay. The significance is defined as \/— 2 ln(£o/'^max), where 



•-max 



and Co are likelihoods for the best fit and zero signal yields, respectively. The significance 
is corrected for systematics due to the peaking background, which is estimated using AE 
and Mbc sidebands. Efficiencies include all intermediate resonance branching fractions 2^ 



and were obtained from MC simulation, assuming = 0"^ for the D*j{2317) and = l"*" 
for the Dsj(2460). We assume equal production of neutral and charged B mesons. 

The PDG value of B{D^ 07r+) = (3.6 ± 0.9)% with a 25% uncertainty is 
used to obtain the branching fractions listed in Table 1. BaBar has recently deter- 
mined the branching fraction B{D^ </)7r+) = (4.81 ± 0.52 ± 0.38)% [i^, which has a 
smaller uncertainty of 13%. If we use this BaBar value the product branching fractions 
become B{W D*j{2317)+K-) x B{D*j{2317)+ D+n^) = (3.3 ± 0.6 ± 0.7) x 10"^ and 
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TABLE I: Signal yields, efficiencies, product branching fractions, and significances for the B — >■ 
DfjK~ processes. The first error is the statistical uncertainty, the second is the systematic un- 
certainty, and the third error is the uncertainty due to Df decay branching fractions. Product 
branching fractions are obtained from simultaneous AM{Dsj) fits of three Dg decay modes as 
described in the text. 



Decay mode 


Yield 


Yield 


Efficiency 


Product B{B^ D^jK-)x 


Signif. 




AM{Dsj) 


AE 


(10-4) 


B{Dsj ^ DsTT^i^)) (10-5) 


a 


D:ji2317)+K- 


35.3 ±6.4 


34.1 ±6.6 


21.9 ±0.6 


4.4 ± 0.8 ± 0.6 ± 1.1 


9.2 


Dsj{2460)+K- 


11.2 ±5.4 


10.2 ±5.4 


59.5 ± 1.4 


0.53 ± 0.20i[!;}^ 
< 0.86(90%C.L.) 


3.1 



B{B^ Dsj{2A60)+K-) x B{Dsj{2A60)+ D+7) = (0.40 ± O.lbtoH) x 10-^. The major 
sources contributing to the systematic error are shown in Table 2. More details about the 
systematic uncertainties can be found in 0. 



TABLE IL Systematic uncertainties in the B — s- D'JjK branching fraction measurements. 





Systematic error (%) 


Source 


Dij{2mrK- 


Dsj{2m)+K- 


Track reconstruction 


±1 X N tracks 


±1 X N tracks 


Charged particle identification 


±2 X Nparticles 


±2 X Nparticles 


Photon and vr*^ reconstruction 


±5 


±2 


reconstruction 


±3 


±3 


AE and likelihood ratio shapes 


±4 


±4 


Helicity angular distribution assumption 


±4 


+9 
-0 


Background subtraction 


±6 


±5 


Fitting procedure 


±3 


±5 


MC statistics 


±2 


±2 


Number of BB pairs 


±1.5 


±1.5 


Total 


±14 


+16 
-13 



In conclusion, improved measurements of — > DfjK~ decay modes have been per- 
formed using a data sample approximately 2.5 times larger. Good agreement with the previ- 
ous measurement is obtained |7|. The value of B{B'^ D*j{2317)^K~) is of the same order 
of magnitude as B{B^ D^K~) and significantly larger than the B^ Dsj{24:60)~^K~ 
branching fraction. 
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